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A NOTE ON KINETIC ENERGY, DISSIPATION AND ENSTROPHY 

JIE-ZHI WU* * * § , YE ZHOU* , AND MENG FAN§ 


Abstract. The dissipation rate of a Newtonian fluid with constant shear viscosity can be shown to 
include three constituents: dilatation, vorticity, and surface strain. The last one is found to make no 
contributions to the change of kinetic energy. These dissipation constituents are used to identify typical 
compact turbulent flow structures at high Reynolds numbers. The incompressible version of the simplified 
kinetic-energy equation is then cast to a novel form, which is free from the work rate done by surface stresses 
but in which the full dissipation reenters. 
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1. Introduction. We make a general theoretical examination on the relation between kinetic energy 
transport, dissipation, and enstrophy. First, we derive an exact but simplified transport equation for the 
kinetic energy of viscous compressible flow, in which the full dissipation function is replaced by the squares 
of vorticity and dilatation. This forms a local counterpart of the classic Bobylcff-Forsyth formula [1] on 
the integral equivalence of incompressible dissipation and enstrophy under special boundary condition, and 
provides a physical clarification on dynamic processes really involved in the evolution of kinetic energy. The 
dissipation constituents are used to identify typical compact flow structures. Second, we cast the incom- 
pressible version of the simplified kinetic-energy equation to a novel form, of which the local spatical average 
docs not explicitly depend on boundary conditions, but the full dissipation reappears along with enstrophy. 
This result is of relevance to some current studies of intcrmittency and scaling laws in inhomogeneous and 
anisotropic turbulence. 

2. Dissipation constituents. Let u = V x u and = V u be the vorticity and dilatation, S and T be 
the strain-rate and stress tensor, respectively, so that for Newtonian fluid T = (— p + \$)I + 2//S (c.g. [1,2]), 
where I is the unit tensor, and ^ and A are the shear and second viscosity 1 . Throughout this note we assume 
fi is constant. The conventional transport equation for kinetic energy E = |ti| 2 /2 reads 

DE 

(1) p~=r- =p/-ti + tfp + V-(T •«)-$, 
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1 One may also replace S by its derivatoric part S' = S — so that T = (— p T Ctf)I + 2jxS', where C = A + |/r is the bulk 
viscosity coefficient. This does not affect our argument but is less convenient. 
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where / is a body force and 4> = At? 2 + 2/iS : S is the dissipation rate. It is straightforward to derive an 
invariant decomposition <£ = $ + , where 

$ = (A + 2p)t? 2 + p\v 2 , 


(2) $ D = -V.(2/iB-tl), 

are, respectively, a positive-definite part caused by t? 2 and |o;| 2 , and a non positive- definite “dissipation” 
caused by surface- strain rate tensor 2 [3] 

(3) B = t?I - Vu r with V ■ B = 0. 

From (2) the compressible extension of the Bo by leff- Forsyth formula follows: 

(4) [ $dV = f [(A + 2p)t? 2 + p\u>\ 2 }dV - 2fi f n B udS. 

Jv Jv JdV 

Thus, if n • (B • u) = 0 on dV, we have an integral equivalence of and 

A similar decomposition of T has been introduced by Wu and Wu [4-5]: 

(5) T = T + T b , T = -IH + 2pf), Ts = — 2/xB, 

where II = p — (A 4- 2/x)t? is the isotropic part of T and fl is the vorticity tensor. Substituting (5) into the 
Cauchy motion equation reveals that T# plays no role in momentum balance. Now, substituting (5) into 
V * (T ■ u) in (1) leads to a splitting of the work rate W per unit volume done by surface stress: 


W = V • (T • u) = W p 4- Wt + W» + W B , 


where 


W p = - V • (pti), W* = V [{a + 2p)t?u], 


W u — - V • (fiu) x «), W# = -V • (2pB * ti). 

Thus, VF/j and <£/? are always self-balanced, making no contribution to the change of E. Hence, Eq. (1) can 
be simplified to 

(6) p~ = pf ■ u + tip + V • (T u) - 

which proves a local “equivalence” of and <£ under the conditi )n that T be replaced by T simultaneously. 
Whenever Tb does a work (positive or negative), it is directly and locally dissipated to <3 >b which, as is 
easily proved, does cause a change of the internal energy. Therefore, so far as its mechanical aspects arc 
concerned, a Newtonian fluid can well be simplified to a hypothetic medium of which the motion consists of 
only an isotropic expansion/compression and a spin. A significant saving of this simplification in numerical 
computation was cited in [5]. The surface deformation, which is generically most complicated, plays no role 
except on free surface [4]. 

2 In this report and Dishington’s paper there is a sign difference in the definition of B. 
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In the above analysis a key issue is the distinction between and $>b or the relation between 4> and 
<I>. Studying their distribution characters allows us to identify different coherent and compact structures in 
high Reynolds-numbcr flows. Typically, these include two shearing structures: thin vortex layers and tubes, 
and a compressing structure, shock waves. While across a shock there is evidently an extremely strong peak 
of tf 2 , the distribution character of c o 1 and for shearing structures can be easily illustrated by idealized 
models. For a unidirectional shear flow with u = (U (y, t), 0, 0) in Cartesian coordinates (x, y , 2), there is 

$ = $ = nu) 2 = ( , $ 0 = 0 ; 

while for a pure vortex with u — (0, v(r, £), 0) and u = (0, 0,o;(r, £)) in cylindrical coordinates (r,0, 2), there 
is 

^ 9 i; dv 

4? = (MjJ , = W B = — 4/x- — , $ = j-i 

r or 

Note that inside the vortex core (r < r 0 ) with dt'/dr > 0 there is a negative peak of which cancels the 
peak of (<£> = $ ma x at r = 0), so that 4> remains small as it should (4> = 0 in a solid-like core). But 
outside the core reaches a hollow-tube like (macaroni- like) positive peak $ B max at an r = r m > ro, 
which becomes the main contributor to <h. For r >r m the contribution of enstrophy to $ is ncglegible. The 
peak behaviors of $ and for layers and tubes can be easily checked by considering, say, U (y) = tanh(i//^) 
with S <C 1 and an Oseen vortex with //<1, respectively. 

When Re^> 1, these idealized models can serve as local building blocks of more generic layer and tube 
structures, provided that the strain components along the layers or tube axes are weaker than the vorticity 
therein. In particular, in a compressible turbulence the strong peaks of a; 2 , and the simultaneous peaks of 
u) 2 and $ arc good indicators of random shocklets, filaments, and shearing sheets, respectively. Indeed, in a 
numerical simulation of two-dimensional supersonic turbulence, Porter et al. [6] have observed complex sheet- 
like shocklet structures by visualizing dilatation field. From their direct numerical simulations of a uniformly 
sheared incompressible turbulence, Tanaka and Kida [7] used the ratio 3>/4> to identify many thin vortex 
filaments with $/<£> > 2 (as shown above, these high-enstrophy structures must be surrounded by weaker 
hollow-tube peaks of strain), while at the high enstrophy region with comparable strain 1/2 < <$/<!> < 4/3, the 
vortex sheet structures were observed. These results are confirmed by Boratav and Pclz [8], who conducted 
direct numerical simulations of an unforced flow and investigated the strain-vorticity correlation in a plane 
spanned by / fj, and enstrophy 4>///. Now, our result provides a theoretical support for these identifications. 
Interestingly, for the Oseen vortex at t — 1 we found ® B max/$max = 0.Q9, completely independent of the 
Reynolds number. The ratio r m /ro, which is about 1.2, has only a very weak dependence on Re. Thus, 
as Re — ► 00 and r 0 — > 0, the distribution of $ and (or the full 4>) for a filament becomes a singular 
“spaghetti” enclosed by a singular “macaroni”, having the same fractal dimension (approaching 1). 

3. Kinetic-energy equation. We now further focus on incompressible flow with p = 1, which is a 
typical situation where the thermodynamics is not involved. In this case (6) reduces to 

(7) ^ = / • u - V • (pu) - j/V • (w x ti) - $. 

In the studies of turbulent-energy budget one removes the divergence terms in (7) by taking average over a 
region with periodic boundary conditions. Then for statistically steady turbulence and as v — > 0, one obtains 
a balance solely between the work of stiring force / and energy cascade in integral range, and a balance solely 
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between the energy cascade and dissipation (enstrophy) in ineitial range (e.g. Frisch [9]). This procedure 
is necessary to obtain the Kolmogorov scaling law and its mod< rn revisions. However, for the local scaling 
of structure functions in an inhomogeneous and anisotropic turbulence, the boundary condition on a small 
sphere (in which the average is to be taken) must be kept generic. Therefore, it will be highly interesting 
if a general incompressible energy equation can be cast to a form free from work- rate terms, so that the 
averaging process does not explicitly rely on specific boundary conditions. We now consider this issue. 

The advection of E and pressure work in (7) can be combined to V • (ti/io), where ho is the stagnation en- 
thalpy. To remove this term we split the momentum equation inio two parts: a transverse or divergence- free 
part (including harmonic components), and a longitudinal or eurl-frec part. Namely, we make a decom- 
position of / and the Lamb vector / = u x ti, denoted by / == f ± + f\\ and l = l± + Z|j. It then follows 
that 

du 

(8) = f ± -l ± - i/V x u, Vh 0 = -I,,, 

where l± is the projection of advection ti • Vu onto the solenoidal space. From (8a) there is 

BE 

(9) = 

Compared to (7), the pressure work is absent, and advection retains a “residue” l± • u, which is implicitly 
boundary-dependent since the projection is a global operator. Then, since 

$ = v\u)\ 2 , = 2z/V • (u • Vti) = — 2^V 2 p, 


there is 

-i/V • l = zA7 2 /i 0 = iA7 2 £ - ($ - $). 

Therefore, from (9) a novel diffusion equation for the kinetic energy follows: 

(10) ^-t/V 2 )E=(/ x -/x) -«-!(* + $), 

in which dissipation and enstrophy appear simultaneously and symmetrically, but all boundary-dependent 
work-rate terms disappear. Note that the reappearance of $ in ( 10) docs not conflict with (7); it comes from 
a part of the vorticity work (the other part being the diffusion o: E). Moreover, compared to the usual form 
of the energy cascade [9], we now see that the cascade is neatly icpresented by the statistic average of l± ■ ti, 
which is dominated by the transverse Lamb vector and has a compact support where l± ^ 0. 

4. Conclusion. It has been noticed that the relative contributions to the dissipation from strain- 
dominated and enstrophy-dominated structures may be relevant to the strength of intcrmittency, as well 
as to the scalings in longitudinal and transversal structure fur ctions [10]. Now, by placing the enstropy 
and dissipation in equal footing and exhibiting their interplay, Eq. (10) clearly emphasizes that, if a more 
boundary-independent (and hence more universal) theory is to te reached, both strain rate (via dissipation) 
and vorticity (via enstrophy) are important in the energy trai sfer process. On the other hand, if (7) is 
used, one need reconsider the respective roles of full dissipation and enstrophy. In either case a more careful 
treatment of the effect of boundary condition on the average over a small sphere in inhomogeneous and 
anisotropic turbulence is necessary. 
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